Charge carrier induced lattice strain and stress effects on As activation in Si 
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We studied lattice expansion coefficient due to As using density functional theory with particular 
attention to separating the impact of electrons and ions. Based on As deactivation mechanism under 
equilibrium conditions, the effect of stress on As activation is predicted. We find that biaxial stress 
results in minimal impact on As activation, which is consistent with experimental observations by 
Sugii et al. [J. Appl. Phys. 96, 261 (2004)] and Bennett et al. [J. Vac. Sci. Tech. B 26, 391 
(2008)]. 



Stress effects are of great interest in modern ULSI tech- 
nology since they can be employed to improve various 
material properties. Uniaxial stress has been employed 
in MOSFET devices since the 90 nm node technology 
step to improve carrier mobility 1]. Properly applied 
stress can also suppress dopant diffusion [2|, y, |4|, en- 
hance activation [3, la, la| > ancl reduce the band gap |7| . 
Therefore, understanding stress effects becomes essential 
for further MOSFET scaling. 

As deactivation is governed by As m V„ cluster forma- 
tion, and clusters with m = 1 — 4 and n = 1 are con- 
sidered as the dominant species in deactivation kinet- 
ics [8]. Under equilibrium conditions, the concentrations 
of defect X (As, V or As m V) is determined by the free 
As and V concentrations and cluster formation energies: 
Cx = Aexp I —E x /kTj, where E x is the formation en- 
thalpy, and A includes the configuration and formation 
entropy The total chemical As concentration is given by 



C total s~i 
As — ^As 



4 



juCas 



(1) 



Table U lists the formation energies of As m V complexes 
based on the total free energy of 64-atom (or 63-atom, 
with vacancy) super-cells using the density functional 
theory (DFT) code VASP @ with PW91 GGA poten- 
tial [lfj(. All calculations were done at a 250 eV en ergy 
cut-off with 2 3 Monkhorst-Pack k-point sampling [ll| . 
Each time an As atom is added to a vacancy, the forma- 
tion energy is lowered by about 1.5 eV, and thus a larger 
complex is more stable than a smaller one. We calcu- 
lated the As rn V concentrations based on the formation 
energies listed in Table |U Since DFT GGA underesti- 
mates the vacancy formation energy by about 1 eV |l2|, 
we also applied a correction for the As m V formation en- 
ergies using experimental values 13j. AS4V has the low- 
est formation energy and becomes the dominant cluster 
under equilibrium conditions. Smaller clusters can be 
formed during epitaxial As-doped Si growth and early 
stages of annealing, and can dominate before full equili- 
bration is reached 1J, |15J , but we restrict our analysis 
to equilibrium conditions. 

Fig. Q] shows the isolated As concentration as a func- 
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FIG. 1: Equilibrium As concentration and AS4V concentra- 
tion as a function of the total chemical As concentration for 
Ef = E c . Solid lines are plotted with correction for vacancy 
formation energy and broken lines are plotted with DFT for- 
mation energies. Smaller clusters don't appear due to low 
concentration. 



tion of the total As concentration. As the number of As 
forming AS4V increases to become a significant fraction 
of free As, the free As concentration starts deviating from 
the total chemical As concentration, which is consistent 
with previous reports [la, llfj . We should note that the 
AS4V formation energy is actually Fermi level dependent 
due to a charge transfer from the Fermi level to the clus- 
ter when AS4V forms. A higher Fermi level results in 
lower cluster formation energies, and thus the AS4V (As) 
curve becomes steeper (flatter) when the Fermi level de- 
pendent formation energy is used. Fig. [1] uses Ep = E c , 
which is appropriate for degenerately doped Si. 

In equilibrium, the change in the concentration of a 
defect X due to stress is given by 



Cx(0) 



exp 



kT 



(2) 



where AE x (e) is the change in the formation energy of 
X due to stress. In the case of the As m V cluster, it is 
given by [l7J ] 



AE f x (e) = -Vo(A? AsmV - mAeAs) ■ C 



(3) 



TABLE I: Formation energy of As m V clusters. When the ex- 
perimental vacancy formation energy is used (4.60 eV [13|). 
formation energies increase by about 1 eV. The experimental 
value of the V formation energy was calculated by subtracting 
the migration barrier (0.26 eV, DFT value) from the activa- 
tion enthalpy (4.86 eV [13] )■ I n the second row, the first value 
is based on the DFT result, and the second is based on the 
experimental V formation energy. 



TABLE IV: Local lattice structure around an As atom in the 
Si lattice compared to atomic spacing in pure Si. 



INN 
2NN 
3NN 



Si 



As" As+ As (cxp)Tl9f 



2.36 
3.86 
4.53 



2.45 
3.87 
4.53 



2.43 
3.86 
4.52 



2.43 
3.87 
4.53 



V AsV As 2 V As 3 V As 4 V 



E f (eV) 



3.59 
4.60 



2.15 
3.16 



0.68 -0.66 
1.69 0.35 



-2.22 
-1.21 



TABLE II: Induced strain for As and As m V complexes. As 
produces small lattice expansion and As m V complexes result 
in lattice contraction. 



As 



V 



AsV As 2 V As 3 V As 4 V 



Ae 0.018 -0.25 -0.21 -0.22 -0.11 



-0.08 



where Vq is the volume of a lattice, AeAs m v (Aca s ) is the 
induced strain due to As m V (As), C is the elastic stiffness 
tensor of Si, and e is applied strain. The induced strain 
can be determined from the energy vs. strain curve. A 
detailed explanation can be found in Ref. [l7|. The re- 
sults are summarized in Table [TTJ We repeated calcula- 
tions using 216 supercell for selected structures and ac- 
quired the same induced strains. 

As shown in Table [HJ As is calculated to give a 
small lattice expansion. However, several authors have 
observed lattice contractions in heavily As-doped Si, 
which the y a ttributed to free electrons in the conduction 
band 18|, [2fJ, |21 1 . In contrast to their conclusion, DFT 
calculations predict a lattice expansion due to free elec- 
trons in the conduction band ( Table [iTl]) . In Cargill et at, 
the total induced strain (Acas — Ptot&iN a s ) is assumed 
to be given by the sum of the induced strain due to ions 
(Ae As + = (3 s i zc N As ) and free electrons (Ae e = l3 e N As )- 
As shown in Table UTTl the calculated induced strain due 
to As has opposite sign to measured value, but the ab- 
solute difference is small and thus its impact on stress 
effects is minimal. However, the reasoning is very dif- 
ferent in each case, which raises a fundamental question 
about the role of electrons: Do electrons cause expansion 
or contraction in the lattice? To answer this question, 
we performed DFT calculations to find equilibrium lat- 



TABLE III: Induced strain due to As, As + , and free electrons 
and holes. The numbers in parenthesis are extracted from 
Cargill et al. [lg]. Note that in spite of longer As-Si bond 
length in SiesAs 4 " supercell (Table HV)| . the lattice undergoes 
contraction. 



As u 



~A^~ 



~w 



Ae 0.018 (-0.019) -0.22 (0.07) 0.22 (-0.09) -0.26 
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FIG. 2: Change in lattice constant due to free charge car- 
riers. The lattice undergoes expansion (contraction) as free 
electrons (holes) are added. Induced strains are obtained by 
finding equilibrium lattice constant of charged supercell with 
various dopants. One electron in a 64-atom supercell corre- 
sponds to 7.8 x 10 20 cm~ 3 . 



tice constants of charged supercells. From the carrier 
concentration vs. change in lattice constant (Fig. [2]), we 
conclude that electrons expand the lattice while holes 
cause lattice contraction. 

The lattice expansion due to electrons raises another 
question about the relation between Si- As bond length 
and the lattice parameter. We looked into the local struc- 
ture around As in Si matrix to answer this question. As 



10 



o 

y 



i ' ' ' ' i 

O— »As 

o- -□ AsV 
O— ©As 2 V 




T=1000C 



_1_ 



_1_ 



i.Ol 



-0.005 




Strain 



0.005 



0.01 



FIG. 3: Stress effects on As and As m V cluster concentration 
under biaxial stress. Note that the two dominant complexes, 
As and AS4V, have minimal stress effects. 
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in accordance with previous experiments 241 |25|. 

In conclusion, by performing DFT calculations of the 
local structure around As in the silicon lattice, we found 
that lattice expansion due to the larger size of an As 
atom is limited to within 3NN distances. The lattice 
contraction in highly As-doped Si can be explained by 
As m V cluster formation rather than free electron as pre- 
viously suggested [18j- The small induced strain due to 
both isolated As and the dominant deactivated cluster 
AS4V results in negligible stress effects on the carrier 
concentration, in accordance with experimental observa- 
tions H El 



FIG. 4: Stress effects on As and As m V cluster concentration 
under biaxial stress. Note that two dominant complexes, As 
and AS4V, have minimal stress effects. 



listed in Table IIVI DFT calculations agree with experi- 
mental measurement up to the 3NN distance and predict 
a local volume expansion around As la, LUJ, |22| . How- 
ever, this expansion is attenuated as distance increases 
and As-Si 3NN spacing is very similar to Si-Si 3NN dis- 
tance. Therefore, changes in the INN bond length are 
not directly linked to changes in the lattice parameter, 
and care should be taken when linking short range atomic 
spacing to lattice constant. In fact, As + produces a lat- 
tice contraction (Ae = —0.22) in spite of longer As-Si 
bond length. A free electron in the conduction/impurity 
band overcompensates this contraction, and thus neutral 
As results in an overall tiny expansion (Ae = 0.018). 

Based on the our analysis, it is likely that experimen- 
tally observed lattice contractions originate from reasons 
other than free electrons. We attribute them to high 
concentrations of vacancies in the form of As m V n clus- 
ters, and find that a vacancy concentration of about 15% 
of the As concentration can reproduce the lattice con- 
traction observed by Cargill et al. 



18] 



This level of 

vacancy concentration was reported based on ab-initio 
calculations by Berding et al. [8| and positron annihila- 
tion spectroscopy by Borot et al. [23j. 

Effects of stress on As and As m V concentrations are 
plotted in Fig. [3] based on Eqs. [2] and [3] The con- 
centrations of the two dominant configurations, As and 
AS4V, undergo changes in opposite directions under biax- 
ial stress, but the magnitude is minimal due to the small 
induced strain. Finally, the free As concentration as a 
function of the total As concentration is plotted in Fig. 2] 
At a given total As concentration, compressive biaxial 
stress enhances As m V formation, and thus the number of 
active As decreases. However, stress effects are minimal 
due to the small induced strains of dominant structures, 
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